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ABSTRACT 

While iron emission lines are well studied in black hole systems, both in X-ray binaries and Active Galactic 
Nuclei, there has been less of a focus on these lines in neutron star low-mass X-ray binaries (LMXBs). How- 
ever, recent observations with Suzaku and XMM-Newton have revealed broad asymmetric iron line profiles in 4 
neutron star LMXBs, confirming an inner disk origin for these lines in neutron star systems. Here, we present 
a search for iron lines in 6 neutron star LMXBs. For each object we have simultaneous Chandra and RXTE 
observations at 2 separate epochs, allowing for both a high resolution spectrum, as well as broadband spectral 
coverage. Out of the six objects in the survey, we only find significant iron lines in two of the objects, GX 17+2 
and GX 349+2. However, we cannot rule out that there are weak, broad lines present in the other sources. 
The equivalent width of the line in GX 11+2 is consistent between the 2 epochs, while in GX 349+2 the line 
equivalent width increases by a factor of ^-^3 between epochs as the source flux decreases by a factor of L3. 
This suggests that the disk is highly ionized, and the line is dominated by recombination emission. We find 
that there appears to be no specific locations in the long-term hardness-intensity diagrams where iron emission 
lines are formed, though more sources and further observations are required. 
Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

Accretion disks around compact objects in AGN and X- 
ray binaries exist in strong gravity around these objects, and 
thus can provide a means to study gravity in the strong- 
field regime. In many of these objects, an Fe K-shell flu- 
orescence emission line is observed in the X-ray spectra 
(e.g. IW hite et al. 1986; Nandraetal. 1997; Revnolds 1997t 
lAsaiet al.ll2000t IMiUer et al.ll2004i) . These lines are formed 
in the inner accretion disk and are observed to be broadened 
due to stro ng Doppler sh ifts and gravitational redshifts in that 
region (see lMilleij|20"07l for a review of iron lines). As these 
lines are shaped by the motions of gas in the disk, they are 
an important diagnostic of the inner accretion disk radius in 
these objects. 

Observations of black holes (in both AGN and X-ray bina- 
ries) hav e revealed so me of these lines to be highly asy mmet- 
ric (e .g. iT^naka et a l. 1995; Fabian etal. 2002; Mille r et al.l 
|2004') due to the extreme motions of gas in the inner disk 
around black holes. While iron lines in black holes are now 
well-studied, iron lines in neutron star low-mass X-ray bina- 
ries (LMXBs) have not been studied in the same detail. Nev- 
ertheless, iron lines in neutron star LMXBs are well known, 
and these lines a re significantly weaker th an in black holes 
(e.g. [White et al.llT985. 19 86; Hi rano et"anil987l: lAsai et all 
l2OO0t iDi Salvo et al.ll2005l) . Thus, flie line shapes could not 



Electronic address: ecackett@umich.edu 

' Department of Astronomy, University of Micliigan, 500 Churcli St, Ann 
Arbor. MI 48109-1042, USA 

- Dean McLaughlin Postdoctoral Fellow 

^ MIT Kavli Institute for Astrophysics and Space Research, MIT, 70 Vas- 
sar Street, Cambridge, MA 02139-4307, USA 

Astronomical Institute 'Anton Pannekoek', University of Amsterdam, 
Kruislaan 403, 1098 SJ, Amsterdam, the Netherlands 

' Kapteyn Astronomical Institute, University of Groningen, P.O. Box 800. 
9700 AV Groningen, the Netherlands 

* Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138, USA 

' Department of Physics, University of Warwick, CV4 7AL, UK 



be studied in detail, and could be modeled as just a Gaussian. 
Only with recent sensitive XMM-Newton and Suzaku obser- 
vations have broad, asymmetric iron lines been observed in 
neutron stars, and these observations highlight the potential 
to use iron lines in neutron stars as a probe of the inner disk 
radius; the inner disk radius is an upper limit on the neutron 
star rad ius (iBhattacharyya & Strohmayerll2007t ICackett et al.l 
l2008ah . 

The use of these lines as a disk diagnostic also allows for 
a test of the standard neutron star accretion picture described 
by color-color diagrams and kilo-Hertz quasi-periodic oscil- 
lations (kHz QPOs). These kHz QPOs are seen to vary in 
frequency as the neutron star changes state and moves around 
the color-color diagram (see van der Klis 2006, for a com- 
prehensive review). The kHz QPOs have frequencies that 
may be associated with the orbital frequency of the inner disk 
and there are several proposed mechanisms for their produc- 
tion ( e g. 'Miller et al.'1998'; 'St eUa et al.|[T999t iLamb & Milled 
1200 U [Abramowicz et al. 2003). By combining color-color 
diagrams, QPOs and iron lines, it may be possible to break 
the degeneracies in the models. For instance, if they are due 
to orbital motions in the disk, then a change kHz QPO fre- 
quency may correspond to a change in the inner accretion 
disk radius. However, in disk resonance models this would 
not be expected. If the inner disk radius can be measured 
independently using iron lines, there is the potential to test 
these models. In fact, recent observations of iron lines in 3 
systems (Ser X-1, GX 349H-2 and 4U 1820-30) show a good 
agreement between the measured inner disk radius from the 
iron lines and the inferred inner disk radius from kHz QPOs 
(when the highest observed kHz QPO frequency is used; 
ICackett et al.ll2008ah . It is therefore important to identify as 
many neutron stars with iron lines as possible, to determine 
whether such a method is possible. The location on hardness- 
intensity (or alternatively color-color) diagrams, where iron 
lines are present, may also be key, and in particular how this 
corresponds to where kHz QPOs are seen. 
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Moreover, iron lines in neutron star LMXBs can pro- 
vide a test of iron line models for black holes. In sev- 
eral objects extremely broad lines are observed from black 
holes. It has been proposed that the extreme width of 
the Unes r equires th at the black holes are spinning (e.g . 
Iwasawa et al. 199& l\yikns et alJ 1200 it iFabian et al.ll2002l: 
Miller et all2002il2004l) . and thus allowing the disk to extend 
closer in than the inner-most stable orbit for a non-spinning 
black hole. Of course, in neutron stars the solid surface of 
the star prevents the accretion disk from extending any fur- 
ther in than that. Any iron lines in neutron star LMXBs 
should therefore be significantly less broad than those in spin- 
ning black holes. Results from recent sensitive observations 
do find that the neutron star lines are narrower than in the 
proposed spinning black holes, supporting that hypothesis 
^Bhattacharyya & Strohmayer 2007; Cackett et al. 2008a). 

As part of the Ch andra HETGS Z/Ato ll Spectroscopic Sur- 
vey (CHAZSS, see lCackett e"tani2008bl for related work) we 
have obtained observations of 6 neutron star LMXBs. Each 
source in this survey was observed at 2 separate epochs using 
the Chandra High Energy Transmission Grating Spectrometer 
(HETGS). Simultaneously, we observed these sources with 
the Rossi X-ray Timing Explorer (RXTE). In this paper we 
present the search for Fe Ka emission lines in these sources 
and we examine how the presence, or lack of iron lines corre- 
sponds to position in hardness-intensity diagrams. In Section 
|2]we describe the observations and data reduction. Section [3] 
details our analysis and results, and in Section |4] we discuss 
the implications of the iron lines that we find. 

2. OBSERVATIONS AND DATA REDUCTION 

Persistent neutron star LMXBs are frequently referred to as 
being either atoll or Z sources due to the shape of the tracks 
they tr ace on X-ray color-color and h ardness-intensity dia- 
grams dHasinger & van der Klisl[T989l) . The differences be- 
tween the two types of sources is thought to be due to different 
mass accretion rates, with the Z sources accretin g at higher 
fractions of the Eddington rate than atolls (see van der Klis 
I2OO6, for a review). This picture was recently strengthened 
by observations of th e first transient Z source XTE J 170 1-462 
dHoman et alj|2007l) whose behavior became more like atolls 
as it decreased in luminosity. 

Six neutron star LMXBs are part of CHAZSS. Of 
these sources three are classified as atolls (4U 1636-53, 
4U 1735-44, 4U 1820-30) and three ai-e classified as Z 
sources (GX 17h-2, GX 340h-0 and GX 349H-2). We observed 
each of these objects twice with Chandra and simultaneously 
with RXTE. Details of these observations are given in Table 
[T] The Chandra and RXTE data reduction is described in the 
following sections. 

2. 1 . Chandra data reduction 

During each observation the HETGS spectrum was read- 
out with the ACIS-S array operating in continuous-clocking 
(CC) mode. Even though these sources are bright, the 3 msec 
frame time in CC mode means that in the dispersed spectrum 
the count rate per pixel is not high enough to cause pile-up. A 
gray filter 100 columns wide was used over the zeroth order. 
This prevents telemetry saturation yet still leaves the position 
of the zeroth order visible, which is vital for accurately ex- 
tracting spectra. 

We performed the data reduction for each observation in the 
same manner using CIAO version 3.3.0.1 and following the 



standard analysis threads'^. The data was reprocessed from 
the level 1 event file to use the latest calibrations. The lo- 
cation of the zeroth order was determined using the tgdetect 
tool. Diffracted grating events were assigned to spectral or- 
ders (using tg_resolve_events), with the spectral orders sep- 
arated using the intrinsic energy resolution of the detectors. 
The data was then filtered for bad grades and for a clean sta- 
tus column (status=0) before applying the good time interval 
filters. Finally, before extracting the spectra the de streak tool 
was used to remove streaking effects on the ACIS-S4 chip. 

The tgextract tool was used to extract the +1 and -1 order 
spectra for the High Energy Grating (HEG), with the nomi- 
nal instrument resolution. We used the standard redistribution 
matrix files (RMFs) from the CALDB, and generated the an- 
cillary response files (ARFs) using the fullgarf script. The 
HEG +1 and -1 orders were examined to ensure that there 
was no wavelength shift between the orders before the spec- 
tra and ARFs were combined using the add_grating_spectra 
script, giving a combined first-order spectrum for the HEG. 
We do not use the Medium Energy Grating (MEG) data here 
as the effective area drops off significantly over the Fe K band 
(6.4-6.97 keV). 

2.2. RXTE data reduction 

The RXTE observations we present here are from pointed 
observations using the proportional counter array (PCA). The 
PCA data were reduced using the LHEASOFT suite (ver- 
sion 6.1). We performed the standard 'goodtime' filtering 
and 'deadtime' corrections using data when the pointing off- 
set < 0.02° and earth-limb elevation angle was larger than 
10°. We use data from PCU 2 only as this is the most reli- 
able of all five PCUs. Spectra were extracted from the 'Stan- 
dard 2' mode data, with systematic errors of 0.6% applied 
to each channel. We used the appropriate background model 
depending on whether the source was faint or bright (with 
faint defined as a source count rate of < 40 counts/s/PCU). 
Background subtracted lightcurves were also created from 
the 'Standard 2' mode data. X-ray bursts were searched for 
in the lightcurves of all observations but were only seen in 
4U 1636-53. We remove all data between 20 s before and 
200 s after the X-ray bursts. 

In order to study the spectral state of the sources we chose 
to create hardness-intensity diagrams for each source. For this 
purpose, we define the hard color as the ratio of count rates 
in the 9.7- 16.0 keV and 6.4-9.7 keV bands. We extracted 
lightcurves (with 256s binning for the atoll sources and 128s 
binning for the Z sources) in each of these energy bands to 
produce hard color curves that we normalized using simu- 
lated Crab-like spectra, to account for long term changes in 
the PCA response. These Crab-like spectra were produced 
with the "fakeit" command in XSPEC, using response matri- 
ces that were created by PCARSP VI 0.1. We also extracted 
the source intensities across the 3.0-25.0 keV band, which 
are also normalized in the same way. 

As we want to compare each source with its position in 
the long-term hardness-intensity diagram, we obtained every 
RXTE/PCA observation of each source in our survey that was 
publically available in the archive before January 2007. We 
reduced the data in the same manner as above to produce the 
long term hardness-intensity diagrams. 

3. ANALYSIS AND RESULTS 

http://cxc.harvai'd.edu/ciao/threads 
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3.1. Spectral Fitting 

In observing iron lines it is important to have a well-defined 
continuum on either side of the Fe K band (6.4-6.97 keV). 
We therefore fitted both the Chandra HEG and RXTE PCA 
spectra simultaneously. The addition of the RXTE PCA spec- 
tra provides a larger energy range above the Fe K band where 
the continuum is well-defined (as opposed to fitting the Chan- 
dra data alone), and thus allows for a more accurate deter- 
mination of the model components necessary to describe the 
continuum emission, facilitating the search for iron lines. 

For a long time, there has been ambiguity as to the cor- 
rect spectral model to fit to neutron star LMXBs, with a 
variet y of different spectral models all fitting equally well 
(e.g.. White et al.|[T98 8: Mitsuda et al. 1989). Recent work by 
iLin. Remillard. & Ho man (2007) explored the different mod- 
els by studying how the measured temperature of the thermal 
components varied over a wide range in luminosity. This was 
done using extensive RXTE observations of two transient atoll 
sources Aql X-1 and 4U 1608-52. Their preferred spectral 
model was determined from the model where the measured 
temperature best followed L (x T*, and differs for the hard 
state and the soft and transitional states. For the hard state, 
they find that a single temperature blackbody plus a broken 
power-law performs best, with the broken power-law used to 
mimic weak Comptonization. For the soft and transitional 
states, their preferred model has two thermal components (a 
blackbody and a multi-color disk blackbody) representing the 
boundary layer and accretion disk, plus a constrained broken 
power-law to model Comptonization. The constrained broken 
power-law has the break energy fixed at 20 keV, and the first 
power-law index constrained to be less than 2.5. 

3.1.1. Hard state observation 

Here, we adopt a very similar prescription to iLin et alj 
(^007). Of all the observations, we find that only one is in 
the hard state (the first observation of 4U 1636-53). Sp ectral 
fitting was performed using XSPEC vll (lArnaudlll996i) . For 
this hard state observation we modeled the continuum using 
a single-temperature blackbody (bbody in XSPEC) plus a 
broken power-law (bknpower in XSPEC), absorption from 
the interstellar medium was accounted for using the model 
phabs. We fitted the Chandra HEG spectra in the 2-9 keV 
range and the RXTE PCA spectra between 3.6 and 25 keV. In 
the RXTE/PCA spectrum above 25 keV the background dom- 
inates over the source flux, we therefore do not fit any ener- 
gies higher than this. Below 2 keV we find that the Chandra 
HEG spectrum turns up slightly, regardless of the model being 
used. We attribute this to scattered light off the telescope sup- 
port structure which is not quite corrected for with the current 
calibration. 

When performing the spectral fitting, we tie all the con- 
tinuum parameters between the Chandra and RXTE spectra 
except the absorbing column density, which was allowed to 
differ. We found that this significantly improved the fit, help- 
ing to correct for the scattered light. Importantly, this does 
not affect the fit through the Fe K band and at higher ener- 
gies. A constant is also included in the model to allow for dif- 
ferences in absolute flux calibration between the two instru- 
ments. Moreover, we found that a broad iron line was present 
in the RXTE/PCA data, thus we added in a Gaussian compo- 
nent (with a central energy constrained to be within 6.4-6.97 
keV) to account for this. However, it is the properties of any 
iron line present in the Chandra data that is of interest here 



- the significantly higher spectral resolution of the Chandra 
gratings is required to study the iron line properties in de- 
tail. The spectral fitting results are given in Table |2] The 
fluxes given are from fits to the Chandra data, with the 10-20 
keV flux determined from an extrapolation of the best-fitting 
model to higher energies. 

3.1.2. Soft/Transitional state obser-vations 

All other observations are in the soft/transitional states, 
thus we use a blackbody (bbody) plus a disk blackbody 
(di skbb) model. When statistically required we also include 
a power-law component, constrained to have a spectral index 
less than 2.5. Again, we fitted the Chandra HEG spectra in 
the 2-9 keV range, however, for the RXTE PCA spectra, we 
only fit f r om 3. 6 to 20 keV. This allows us to simplify the 
iLin et al.l (l2007h prescription - their model has a constrained 
broken power-law with the break energy fixed at 20 keV. By 
ignoring data above 20 keV we can simplify the model by just 
using a single unbroken power-law without affecting the fit at 
lower energies. For the two brightest sources (GX 17h-2 and 
GX 349H-2) we found that the scattered light effect was worst, 
and hence we had to restrict the energy range of the Chandra 
HEG to 3 - 9 ke V for those spectra only. 

Again, we included a constant to account for differences in 
the absolute flux calibration between the instruments. As with 
the hard state observation, all continuum parameters were tied 
between the Chandra and RXTE spectra except the absorb- 
ing column, which was allowed to vary (for the same reasons 
mentioned above). In some cases where the column density is 
not well-constrained and tends to zero in the Chandra spec- 
tra, we fix the column density at the Galactic value in that 
direction, while still allowing it to be a free parameter in the 
RXTE data. A Gaussian component was used to model the 
broad iron line seen in all RXTE spectra apart from GX 340H-0, 
where one was not required. We stress again that the limited 
spectral resolution of RXTE/PCA hampers the use of that data 
to study the iron line properties in detail. 

We found that for all the atoll sources in our sample 
(4U 1636-53, 4U 1735-44, and 4U 1820-30) the constrained 
power-law component was statistically required. However, 
for the Z sources (GX 17h-2, GX 340H-0 and GX 349H-2) the 
constrained power-law did not significantly improve the qual- 
ity of the fit, and thus was not used. The spectral fitting re- 
sults are given in Table [51 and we show the Chandra HEG 
and RXTE PCA spectra for the first observation of GX 17h-2 
in Fig. [U along with the best-fitting model components to the 
Chandra spectrum. 

3.1.3. Searching for iron emission lines in the Chandra spectra 

Having determined the model components needed to fit the 
continuum, we look for iron emission lines in the Chandra 
HEG spectra. This was performed by fitting a Gaussian to 
the Chandra HEG spectrum with the energy constrained to 
be within the 6.4-6.97 keV Fe K range, and determining 
whether this additional model component improves the fit sig- 
nificantly. Note that after adding in the Gaussian component 
to the model, we still allow the continuum parameters to vary 
in the fit. We detect significant iron lines in the Chandra HEG 
spectra of two objects: GX 17h-2 and GX 349H-2. In these ob- 
jects an iron line is present in both observations. Significant 
iron lines are not observed in the Chandra spectra of the other 
four objects. We estimate the line significance from the un- 
certainty in the equivalent width. The detected iron lines are 
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shown in Fig. |2] and the parameters from the Gaussian fits are 
given in Table |4] 

We put upper Hmits on Hnes in the other spectra by fitting a 
Gaussian with a fixed width of ct = 0.1 keV and also a = 0.5 
keV. We choose these line widths to probe for both narrow 
and broad lines. The energy of the line was allowed to be 
a free parameter, apart from those cases where it could not 
be constrained within 6.4-6.97 keV. In those cases, the line 
energy was fixed to 6.7 keV, consistent with the detected lines. 
The determined upper limits are given in Table |5] 

The iron lines observed in GX 17+2 and GX 349+2 do not 
appear to be asymmetric. This is due to the low signal to noise 
in the continuum redward of the line, as Suzaku observations 
of the Une in GX 349 +2 clearly show it as being asymmetric 
dCackett et al.ll2008al) . A simple Gaussian fits the lines well 
in aU observations. While the more complex r elativistic ac- 
cretio n disk emission line model, diskline jFabian et alj 
fT989h . can also fit the lines well, when all the parameters 
are left free in the fit they are not well constrained and the 
model does not statistically improve the fit over a Gaussian. 
To demonstrate that an asymmetric line cannot be ruled out 
we show both the Chandra and Suzaku data together in Fig. 
[3j The inner disk radius determined from fitting the Chandra 
GX 349+2 data with the diskline model can be signifi- 
cantly improved by using the emissivity index, /3 and disk in- 
clination, /, found from fits to the Suzaku data in Cackett et al. 
( l2008ah . We therefore fitted the 2nd Chandra observation 
(where the iron line is the strongest) with the diskline 
model, fixing [3 = -4.1, / = 23° and constraining the line en- 
ergy to be within 6.4-6.97 keV. Doing this, we measured the 
inner disk radius, 7?,„ = 10.45 ±0.49 GM/c^, normalization 
= (1.1 ±0.4) X 10-2, and equivalent width = 92 ± 31 eV. The 
equivalent width is larger than fitting a Gaussian due to the 
asymmetric profile of the diskline model. 

The equivalent width of the lines observed in GX 349+2 
are smaller than those observed with Suzaku (ICackett et alJ 
|^08a), where they find EW = 76 ± 6 eV. However, this differ- 
ence can be attributed to the asymmetric line profile given that 
fitting the 2nd observation with a diskline model gives a 
larger equivalent width than fitting with a Gaussian, and one 
that is consistent with the Suzaku value. We note the first ob- 
servation is consistent wit h the Gaussian fit to the line in this 
object byl laria et all (I2004') . who find 39 ± 13 eV. The equiv- 
alent widths measured for GX 17+2 are consist ent with previ- 
ous observations of the iron line in this source dFarinelli et al.l 
12005; Di Salvo et al. 2000). 

Comparing the equivalent widths of the detected lines be- 
tween each epoch we find that the lines in GX 17+2 are al- 
most consistent with remaining unchanged between the ob- 
servations. On the other hand in GX 349+2 the equivalent 
width is larger by a factor of ^3 in the second observation, 
though this increase is of low significance. One might expect 
that the line equivalent width would increase with increasing 
hard X-ray flux (with energies greater than the line energy). 
However, we find that the 10-20 keV flux decreases between 
observations, the opposite of the change in equivalent width. 
The 2-10 keV flux also decreases between observations by 
a similar amount, thus the overall spectral shape is not very 
different. 

Of the sources where we do not detect iron lines in the 
Chandra spectra, only 4U 1820-30 has had a previously de- 
tected iron line with CCD-resolution spectroscopy. Recent 
Suzaku observations detected a broad line in 4U 1820-30 



( ICackett et alj|200^ab . Moreover, a Une was also ob served 
in the object previously with ASCA (lAsai et al.ll2000h . Nev- 
ertheless, we do not detect any line here. However, the upper 
limits we determine with a = 0.5 keV are not too far from the 
eq uivalent width as m easured in recent Suzaku observations 
bv lCackettetaP (l2008a) . where EW = 51 ± 1 1 eV. This indi- 
cates that the lower signal to noise here may prevent a signif- 
icant detection of the line. We therefore cannot exclude that 
weak, broad lines may also be present in the other sources as 
the observations are not sensitive enough to detect them. In 
fact, we note that in 4U 1636-53 there is evidence of a broad 
iron line in recent Suzaku observations (Homan et al, private 
communication), a s well as from XMM-Newton observations 
jPandel et al.ll2008h . 

3.2. Hardness-intensity diagrams and timing analysis 

The location of the CHAZSS observations on the long-term 
RXTE hardness-intensity diagrams is shown in Fig. |4] The 
first CHAZSS observation is in red, the second CHAZSS ob- 
servation in blue and the other observations in black. As dis- 
cussed above, significant iron emission lines are observed in 
the Chandra spectra of GX 17+2 and GX 349+2, but not in 
any of the other sources. Comparing the hardness-intensity di- 
agrams it appears that there may be no special location where 
iron emission lines are preferentially observed (though obser- 
vations in more states are needed to support this). The impli- 
cations of this are discussed in the next section. 

Of our observations of the two sources that show iron lines, 
only the first observation of GX 17+2 is in the correct state 
to display kHz QPOs (as indicated by its location on the 
hardness-intensity diagram). To investigate this further we 
extracted fast timing mode data from the 3.6 - 60 keV band 
using all available PCUs (which ranged from between 2 and 4 
during the observations). We created a broadband power den- 
sity spectrum using the standard fast Fourier transform tech- 
niques (vanderKlis 1989, 1995). For the first GX 17+2 ob- 
servation we find that there is a pair of kHz QPOs present 
in the power spectrum, with frequencies of 594 ± 1 1 Hz and 
906 ±15 Hz, and fractional rms amplitudes of 2.4 ± 0.3% and 
3.3 ±0.4% respectively. The quality factor which is defined as 
Q = iy()/2A (where I'o is the centroid frequency and A the half- 
width at half-maximum) is 5.1 ± 1.2 for the lower frequency 
kHz QPO and 3 .5 ± 1 .0 for the upper kHz QPO. We find there 
are no kHz QPOs in the second observation of GX 17+2. We 
also analyzed the fast timing mode data for both GX 349+2 
observations in the same manner and found that no kHz QPOs 
were present, as expected based on the sources position on the 
hardness-intensity diagram during these observations. 

4. DISCUSSION 

We present Chandra and RXTE observations of six neu- 
tron star low-mass X-ray binaries at two separate epochs 
(4U 1636-53, 4U 1735-44, 4U 1820-30, GX 17+2, 
GX 340+0 and GX 349+2). We detect iron emission lines in 
the Chandra HEG spectra of both GX 17+2 and GX 349+2 at 
both epochs. No iron lines are detected in the Chandra spectra 
of the other four sources, although we do not have the sensitiv- 
ity to rule out that weak, broad lines are present. Comparing 
the position of the sources in hardness-intensity diagrams (us- 
ing the RXTE data) we find that from these observations there 
appears to be no region of the diagram where the lines are 
preferentially observed. However, we note that we only sam- 
ple a small number of states here, and more sensitive observa- 
tions are needed to confirm this. Such a finding is interesting 
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though since we expect that the hardness-intensity diagrams 
are related to the mass accretion rate and the accretion disk- 
corona interaction. For example, changes from bright, soft 
(disk dominated) states to fainter, hard (corona dominated) 
states may indicate changes in mass accretion rate and the 
prominence of the disk or the corona. One might expect to 
see changes in the iron lines as a function of the prominence 
of the disk or corona. 

The 3 Z sources observed are seen across most locations 
in the color-color and color-intensity diagrams. GX 17+2 is 
seen at both the normal branch (NB)/horizontal branch (HB) 
transition as well as the NB/flaring branch (FB) transition. 
GX 340H-0 is seen mainly on the NB, with part of the first 
observation on the FB. GX 349H-2 is observed on the FB in 
both observations. Thus, iron lines are detected in the Chan- 
dra spectra on both the FB, and at the NB/HB and NB/FB 
transitions. 

Previous BeppoSAX observa t ions of GX \l+2 and 
GX 349H-2 (iDi Salvo etal.1 120001 1200 Ih support our finding 
that iron lines are seen in no special location on the hardness- 
intensity (or color-color) diagram. During the BeppoSAX ob- 
servations of GX 17h-2 the source was in the horizontal and 
normal branches. When the data was split based on state, an 
iron line with an equivalent width in th e approximate rang e 
30-40 eV was observed in all spectra ( Di Salvo et al.ll2003 ). 
For the BeppoSAX GX. 349+2 observations the source was ob- 
served on the flaring branch , and both non-flaring and flaring 
spectra showed an iron line jDi Salvo et alJl200Tl) . 



We do not see any iron lines in the Chandra spectra in any 
of the atoll sources in our sample. The majority of the ob- 
servations are spread out over the softest part of the color- 
color diagram - the banana branch, although we do observe 
4U 1636-53 in the hard (island) state. This indicates that we 
do not observe strong lines in either spectrally soft or hard 
states in the atolls. However, we cannot rule out that weak, 
broad lines are present in the spectra as has been seen with 
more sensitive observations of 4U 1820-30 (ICackett et al.l 
12008a!) . In fact, while we don't detect iron lines during the ba- 
nana branch here, observations of other neutron star LMXBs, 
such as the ato ll 4U 1705 -44, do show iron lines during the 
banana branch dPi Salvo"e t al. 2005; Piraino et al. 2007), sup- 
porting that iron lines can be seen in a wide variety of spectral 
states. 

One interesting proposition is that of using sensitive obser- 
vations of iron lines to help constrain models for kHz QPOs. 
Additionally one can combine the measurement of both the 
disk velocity (from iron lines) and orbital frequency (from 
kHz QPOs, ass uming that is their orig in) to determine the neu- 
tron star mass (ICackett et al.ll2008ah For this to be possible 
we would need to see both iron lines and kHz QPOs simul- 
taneously. Here we have detected both an iron line and twin 
kHz QPOs in the case of GX 17h-2. However, the observation 
of GX 17H-2 is not sensitive enough to detect any broad red 
wing of the iron line, and thus we cannot get a reliable mea- 
sure of the inner accretion disk radius by modeling the iron 
line. Nevertheless, this is promising for further studies using 
more sensitive iron line observations, such as can be obtained 
with Suzaka. 

There are several possibilities for the source of ionizing flux 
that gives rise to the iron emission lines. It may be irradiation 



by a corona, by the base of a compact jet, or by the bound- 
ary layer. These observations allow us to investigate these 
possibilities. For instance, the hardest states in the hardness- 
intensity diagram are usually interpreted as being the state 
where the corona is most dominant, thus we may expect to 
have stronger iron lines in harder spectral states. Or, if the 
ionization is too high in the hardest states, the iron lines may 
get stronger as the source gets softer, but would be expected 
to weaken or even disappear in the disk dominated states. In 
addition, the irradiation may be due to the boundary layer, or 
the base of a compact jet. If the irradiation is due to the base 
of a compact jet, then we would expect to see iron lines when 
the jet is present, and not when the jet disappears. 

We find that the iron lines detected in GX 17h-2 and 
GX 349H-2 are both strongest when the sources are closest to 
the NB/FB transition (though the changes in equivalent width 
ar e not highly significant). Th is is similar to what was found 
by iDi Salvo et all ( l2000l 1200 lb for these sources. Comparing 
the 10-20 keV flux with the line equivalent widths, we find 
that in both GX 17h-2 and GX 349H-2 the equivalent width 
increases with decreasing 10-20 keV flux. This can be ex- 
plained if the gas in the disk is very highly ionized already 
and thus the emission is mostly due to recombination. In such 
a scenario a drop in hard flux should lower the ionization lead- 
ing to more recombinations and a stronger line. The fact that 
we see the line energies consistent with the He-like 6.7 keV 
iron line is in agreement with this. The absence of strong lines 
in sources with short periods (such as 4U 1820-30) is also 
consistent with such a picture, as the disks should be highly 
ionized in these sources. 

X-ray and radio observations of neutron star LMXBs 
have associated jet formation and location on the color- 
color/ h ardn ess-intensity diag rar ns (e.g., iPenninx et al.l 
T988t iHiei lming et all [T9903bl: iMighari & Fended 120061: 
Migliari et al. 2007). From a detailed study of GX 17h-2, 
Migliari et al. (2007) find that radio emission is not detected 
in the FB and turns on as the source transitions into the 
NB. The radio flux density increases through the NB and is 
strongest during the HB. We have concurrent radio observa- 
tions of the sources in this study with the VLA (Cackett et 
al., in preparation), and thus can look for any correlations 
between radio flux density and the iron lines. For GX 349H-2 
we do not detect the source at 5 GHz in either observation 
(consistent with its position on the FB), yet for GX 17h-2 we 
have flux densities of 1.80 ±0.05 mJy and 0.9 ±0.07 mJy 
at 5 GHz for observations concurrent with observations 1 
and 2, respectively. Thus, in GX 17h-2 we see higher radio 
flux densities at the NB/HB st ate than in the NB/F B state 
(consistent with the picture of iMigliari et al.l l2007h . Note 
though that while there is a significant change in radio flux 
between observations the equivalent width of the iron line in 
GX 17h-2 does not change significantly. In addition, given 
that we do not detect any radio emission from GX 349H-2 
yet we do see an iron line this would tend not to favor a 
compact radio jet as the source of ionizing flux for the line 
emission. However, it is not clear whether a disk corona or 
the boundary layer is the more likely source of irradiation. 

JMM gratefully acknowledges support from Chandra. 
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TABLE 1 

Details of CHAZSS Chandra and RXTE observations 
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Observation 


Chandra 


Chandra 


RXTE 


RXTE 




Date 


ObsID 


Exposure (ks) 
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91152-06-02 
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91152-02-01 


13.4 


GX 17+2 


2006 Aug 19 


6630 


24.1 
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2006 Oct 20 
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GX 349+2 
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12.6 


91152-01-02 
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GX 349+2 


2006 Aug 20 
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Note. — All Chandra observations were performed with the HETGS/ACIS instrument. 
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TABLE 2 

Best fitting continuum spectral parameters for the hard state 



Object Observation 


Nh (CXO) 


Nh (RXTE) 


Blackbody 




Broken Power-law 




2-10 keV flux 


10-20 keV flux 










kT (keV) Norm. (lO"^)" 




(keV) 


Norm.'' 


(IQ-'ergcm-^s-') 


(IQ-'ergcm-^s-') 




4U 1636-53 1 


0.57±0.19 


1.83 ±0.42 


1.23 ±0.05 0.25 ±0.05 


1.70 ±0.05 


17.1 ±1.6 1.97±0.07 


0.24 ±0.03 


1.1±0.1 


0.6±0.1 


0.80 



Note. — The parameters are from joint spectral fits to the Chandra HEG and RXTE PCA spectra. In the joint fitting all continuum parameters were tied except for Nh- All 
uncertainties are quoted at the 90% confidence level. 

^ Normalization of the blackbody component is defined as (L/IO^^ erg s ')/(D/10 kpc)^'' Normalization of the broken power-law is defined as photons keV ' cm"^ s~' at 1 keV 



TABLE 3 

Best fitting continuum spectral parameters for soft/transitional states 



Object 


Obs. 


Nh (CXO) 


Nh (RXTE) 


Blackbody 


Disk Blackbody 


Power-law 


2-10 keV flux 


10-20 keV flux 








(10^2 cm-2) 


(10^2 cm-2) 


kT (keV) 


Norm. (lO-^f 


kT (keV) 


Norm.'' 


r 


Norm.'^ 


(10"' erg cm"- s"') 


(10"' ergcm"2 s"') 




4U 1636-53 


2 


0.39dz0.12 


1.51±0.32 


2.20 ±0.06 


1.3±0.1 


1.21 ±0.04 


50.3 ±6.7 




0.62 ±0.13 


2.7±0.1 


0.5±0.1 


0.86 


4U 1735-44 


1 


0.30 (fixed) 


1.38±0.19 


2.46 ±0.02 


2.1±0.1 


1.52 ±0.02 


24.7 ± 1.4 


2 50+*" 


0.40 ± 0.04 


3.2±0.3 


0.8±0.1 


0.84 




2 


0.30 (fixed) 


1.06±0.18 


2.59 ±0.03 


2.7±0.1 


1.74 ±0.03 


18.4± 1.8 


2 50+" '' 
■^■■^ -0.07 


0.29 ± 0.04 


3.9±0.3 


1.1±0.1 


0.92 


GX 17+2 


1 


3.18±0.03 


4.63 ±0.08 


2.66 ±0.04 


7.3±0.1 


1.57 ±0.01 


152.5 ±5.7 




12.1 ±0.9 


2.9±0.1 


1.02 




2 


3. 12 ±0.05 


3.67±0.14 


2.54 ±0.02 


4.4 ±0.1 


1.55 ±0.03 


152.2 ±6.8 






10.5±0.8 


1.8±0.1 


0.92 


4U 1820-30 


1 


0.15 (fixed) 


1.74 ±0.33 


2.36 ±0.03 


4.4 ±0.2 


1.39 ±0.05 


62.1 ±9.5 


2. 17 ±0.20 


0.21 ±0.10 


5.6 ±0.6 


1.4±0.1 


0.91 




2 


0.15 (fixed) 


1.82±0.27 


2.41 ±0.03 


6.0±0.2 


1.47 ±0.05 


73.9 ±6.4 


2.08±0.19 


0.18 ±0.07 


7.9 ±0.7 


2.0±0.2 


1.24 


GX 340+0 


1 


6.41 ±0.05 


6.25 ±0.26 


2.66 ±0.04 


2.4±0.1 


1.63±0.01 


130.0±4.5 






8.2 ±0.6 


1.3±0.1 


0.90 




2 


6.47 ± 0.05 


6.29 ±0.28 


2.54 ±0.03 


2.9 ±0.3 


1.59±0.01 


139.3±5.2 






8.1 ±0.6 


1.4±0.1 


1.05 


GX 349+2 


1 


1.79 ±0.06 


3.13±0.14 


2.57 ±0.01 


7.2±0.1 


1.86±0.02 


88.4± 1.1 






15.9± 1.3 


3.4±0.3 


0.96 




2 


1.53 ±0.02 


2.59 ±0.30 


2.46 ±0.03 


5.8±0.3 


1.67 ±0.03 


105.1 ±5.4 






12.5± 1.2 


2.3 ±0.2 


0.96 



Note. — The parameters are from joint spectral fits to the Chandra HEG and RXTE PCA spectra. In the joint fitting all continuum parameters were tied except for Nh- 
All uncertainties are quoted at the 90% confidence level. For GX 17+2 and GX 349+2 the fits to the Chandra spectra also include a Gaussian to model the iron emission line. 
Best-fitting parameters for these Gaussian are given in Table|4] 

" Normalization of the blackbody component is defined as (L/IO" erg s"')/(D/10 kpc)^'' Normalization of the disk blackbody is defined as (R,„/km)/(D/10 kpc)^ x cos 61, where 
R,„ is the inner disk radius and 6 the angle of the disk.*^ Normalization of the power-law is defined as photons keV"' cm"^ s"' at 1 keV 
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TABLE 4 

Iron emission line properties from fits to Chandra spectra 



Object 


Observation 


Energy (keV) 


FWHM (keV) 


EW (eV) 


Line flux (10 ' photons cm ^ s ') 


Significance (u) 


GX 17+2 


1 


6.67 ±0.08 


0.31 ±0.12 


14.1 ±5.2 


2.0 ±0.7 


4.5 




2 


6.70 ±0.06 


0.45±0.17 


33.0±10.3 


3.7±1.2 


5.3 


GX 349+2 


1 


6.75 ±0.06 


0.26±0.14 


18.1±7.2 


3.2±1.3 


4.1 




2 


6.71 ±0.06 


0.52±0.17 


53.0±14.4 


6.8±1.8 


6.1 



Note. — All uncertainties are quoted at the 90% confidence level. 



TABLE 5 

Upper limits on iron emission line non-detections from Chandra 
observations 



Object 


Observation 


Assuming cr = 


0.1 keV 


Assuming u 


= 0.5 keV 






Energy (keV) 


EW (eV) 


Energy (keV) 


EW (eV) 


4U 1636-53 


1 


6.7 (fixed) 


< 12.1 


6.7 (fixed) 


< 80.1 




2 


6.7 (fixed) 


<8.1 


6.7 (fixed) 


< 31.9 


4U 1735-44 


1 


6.76±0.11 


< 26.0 


6.7 (fixed) 


< 39.6 




2 


6.7 (fixed) 


< 14.4 


6.7 (fixed) 


< 29.4 


4U 1820-30 


1 


6.7 (fixed) 


< 8.5 


6.7 (fixed) 


< 31.9 




2 


6.7 (fixed) 


< 11.7 


6.7 (fixed) 


<40.6 


GX 340+0 


1 


6.66 ±0.19 


< 11.2 


6.7 (fixed) 


< 16.2 




2 


6.7 (fixed) 


<9.2 


6.7 (fixed) 


< 8.2 



Note. — Upper limits are quoted at the 95% confidence level. 




Energy (keV) Energy (keV) 

Fig. 1. — Left: Chandra (black) and RXTE/PCA (red) spectra from the first observation of GX 17+2. The bottom panel shows the residuals (x = (data- 
model)/(T). Right: The best fitting model components to the Chandra data are shown (extended to show the model up to 20 keV). The solid line is the overall 
model spectrum, the dotted line is the disk blackbody component, the dashed line is the blackbody component and the dash-dot line is the Gaussian modeUng the 
iron line. 
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Energy (keV) Energy (keV) 

Fig. 2. — Iron emission lines detected in the Chandra HEG spectra of GX 17+2 (top) and GX 349+2 (bottom). The first observations are on the right, and the 
second observations are on the left. Plotted is the ratio of the data to the continuum model. The solid fine is the best-fitting continuum + line model. 
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Fig. 3. — A comparison of the iron emission lines detected from GX 349+2 with Chandra (black, 2nd observation) and Suzaku (red, only XIS detector 2 is 
shown). 



Iron lines in neutron star low-mass X-ray binaries 



11 



1.2 
1 



S 0,8 

'SO.6 

si 



1 1 1 I I I I 1 1 I I I I 1 1 I 1 1 1 
rGX 17 + 2 - 



I I I I I M I I I I I I I I I I I I 

GX 340+0 



0.2 Ft ■ f f I I f > I I I I I I I I I ■ TTi I I I I I I I I I I I I I I I I I i+i I I I I I I I 1 I I i 1 f f I I I r 



n M I I M I I I I I t [ I I I I 
GX 349+2 



0.5 1 1.5 



0.2 0.4 0.6 O.fl 



0.5 I 1.5 2 



1.4 

- 1 

o 

S 0-6 
X 



0,4 r 



J 1 1 1 1 r 1 1 1 1 1 1 1 i_ 

: 4U 1B36-53 : 

M 


1 1 1 1 1 1 1 1 1 1 1 1 M 1 1 1 M 1 1 1 1 1 

- 4U 1735-44 I 


1 1 1 1 1 1 1 p 1 1 1 « 1 1 

" 4U 1820-30 


■■^•^^^^^■^^ ^ ' 

"i 1 1 1 1 1 1 1 1 1 1 1 1 r 


T 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 r 1 iT 


~i 1 1 1 1 1 1 1 1 1 1 1 1 1 i~ 



0,1 0.2 
intensity 



0.1 0.2 0,3 0.4 
Intensity 



0.2 0.4 0,6 0,8 
totensity 



Fig. 4. — Long-term RXTE hardness-intensity diagrams for all six sources. The first CHAZSS observations are red and the second CHAZSS observations are 
blue. The intensities are normaUzed to the Crab in the 3.0-25 keV band. 



